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The  Division  of  a  Circle  or  Spherical 

Surface  Into  Equal- Area  Cells 

or  Pixels 


1.  INTRODUCTION 

The  task  is  to  divide  a  circular  area,  such  as  a  PPI-radar  scope  picture,  into  equal  units 
of  area  or  cells.  Although  the  units  could  be  few  and  large,  generally  an  area  must  be 
partitioned  into  many  small  cells  or  pixels.  If  the  pixels  can  be  made  sufficiently  small,  then 
mapping  a  phenomenon  such  as  areal  rainfall  becomes  a  matter  of  assigning  a  single  value 
(for  example,  rainfall  amount)  to  each  pixel. 

A  flat  circle  can  be  used  for  mapping  properties  such  a  clouds  on  a  sky  dome,  or  properties 
on  the  spherical  surface  of  the  earth.  *  Rossow  and  Carder  addressed  this  problem  and 
developed  and  approximate  solution  for  the  equal-area  net  grid.^ 

The  procedure  presented  in  this  report,  however,  yields  exact  equal-area  pixels. 

2.  METHOD 

The  proposed  method  begins  by  dividing  the  circle  (Figure  1)  of  overall  radius  R  into  a 
central  unit  circle  of  radius  Sr  and  M  rings,  each  ring  of  thickness  26r,  where 

(2Af-rl)^r  =  R  (1) 


The  unit  area  SA  is  that  of  the  central  cell 

SA  =  7r(^r)*. 


(2) 


The  circle  that  includes  the  first  ring  has  radius  SSr,  and  the  area  Ai  within  this  circle 
is  given  by  Ai  =  97r(5r)*. 

The  area  in  the  first  ring,  or  annulus,  =  8x(^r)*.  If  the  ring  is  now  divided  into  8 
equal  cells,  each  cell  will  be  equal  in  area  to  the  central  unit  cell  SA. 
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Figure  1:  Division  of  a  Circular  Area  into  Circles  and  Then  Ekjual-Area  Pixels,  Showing 
Numbering  of  Pixels 
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Likewise,  the  second  ring  will  have  outer  radius  5Sr  and  it  can  be  divided  into  2x8  cells, 
each  having  area  SA.  Similarly  the  mth  ring  can  be  divided  into  8m  cells,  each  having  the 
area  SA.  The  total  number  of  cells  /m  up  to  ^nd  including  the  mth  ring  will  be 

/„  =  (2m  +  l)"  (3) 

for  a  total  area  Am  of 

Am  =  (2m  +  l)^7r(6r)^.  (4) 

Figure  1  depicts  the  division  of  a  circle  into  five  rings,  for  a  total  of  121  equal-area  cells. 
The  shapes  of  the  cells  begin  with  the  small  circle  at  the  center:  cell  number  1.  In  the  first 
ring  the  8  cells  are  shaped  like  pieces  of  pie  with  the  tips  eaten  away.  In  the  remaining  rings 
the  shapes  of  the  elementary  units  converge  onto  the  rectangular  shape  of  length  '2Sr  and 
width  Tr{Sr)f2. 

The  cells  are  numbered,  beginning  with  the  number  1  for  the  central  cell,  and  increasing 
by  8m  with  the  mth  ring  (m  =  1,2,3,...). 


Problem  1:  To  find  a  cell  number  I ,  (fiven  the  polar  coordivalc.s  (r,  0),  when  the  men  mental 
length  Sr  is  known. 

The  cell  number  /  is  determined  as  follows: 


Where  the  symbol  INT{x)  denotes  the  integer  part  of  the  number  x  always  less  than  i), 
then  the  number  of  the  last  whole  ring  n  is  given  by 


n  =  INTin') 


(5) 


and  n'  is  defined  by 

(2n'  4-  l)6r  =  r. 

(6) 

The  number  of  the  ring  in  which  (r,0)  is  located  is  (n  4-  1  )• 
the  n  rings  is 

The  total  number  of  cells  in 

/,  =  (2n4-l)*. 

(7) 

The  number  of  cells  le  remaining  in  the  (n  4-  l)th  ring  is 

h  -  iNT{i;)  +  I 

(8) 

where 

Ig  =  {theta /2ir){8n  4-  8) 

(9) 

if  6  is  in  radians,  and  =  (0/36O)(8n  4-  8)  if  5  is  in  degrees. 


Finally  the  cell  number  I  in  which  (r,  0)  is  located  is: 

/  =  /r  +  /S 


(10) 
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Problem  2:  To  find  the  polar  coordinates  (r,d)  of  the  middle  of  a  cell,  given  the  cell  number 
I  and  the  incremented  unit  of  length  Sr. 


The  number  of  whole  rings  n  is: 
where  it'  is  given  by 


n  =  INTin') 
i2n'  +  lf  =  I 


(Note:  If  n'  is  a  whole  number  then  make  n  =  n'  —  1), 


The  radial  distance  r  to  the  central  point  of  the  cell  is 

r  =  (2n  +  2)Sr 

The  number  of  cells  U  up  to,  and  including,  the  nth  ring  is 

Ir  =  (2n  +  1)* 

Hence  the  remaining  number  of  cells  le  is 

h  =  I -Ir 

which  places  the  central  point  of  the  last  cell  at  angle  0  such  that 

<?/27r  =  (/«-  l/2)/(8n  +  8) 

if  0  is  in  radians. 


(11) 

(12) 

(13) 

(14) 

(15) 

(16) 


If  angular  measurements  are  made  in  degrees,  the 

tf/360  =  (/ff  -  l/2)(8n  +  8).  (17) 

Equations  (13),  (16),  or  (17)  provide  the  coordinates  (r,0). 

Problem  3;  Application  to  the  earth’s  northern  hemisphere. 

This  procedure  can  be  applied  readily  to  a  spherical  surface  which  can  be  projected  onto 
a  flat  circular  equal-area  map.  For  the  northern  hemisphere  the  map’s  center  represents 
the  north  pole  and  the  outer  perimeter  represents  the  equator.  In  the  Lambert  Azimuthal 
Equal-area  projection  the  parallels  of  latitude  appear  as  concentric  circles  centered  at  the 
north  pole,  and  the  meridians  of  longitude  appear  as  straight  lines  radiating  outward  from 
the  center. 


Let  be  the  radius  of  the  spherical  earth.  That  is, 

Hr  =  3437.75nm  =  6371.00km. 
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The  area  of  the  northern  hemisphere  A  is 


A  =  2TrRl.  (18) 

The  equal-area  map  (given  a  scale  factor  of  1.0)  must  have  an  overall  radius  R  such  that 

R  =  ReV2  (19) 

If  we  choose  to  divide  the  circular  map  into  M  rings,  the  incremental  value  8r  is  defined  by 

(2M -\-\)8r  =  R.  (20) 

The  mth  ring  will  have  and  outer  radius  such  that 

r,„  =  (2M  -t-  l)6r.  (  21 ) 

The  map  radius  r  that  corresponds  to  a  parallel  of  latitude  (f>  is: 

r  =  R,^2{\  -  sin  4>).  (22) 

The  map  angular  coordinate  0  that  corresponds  to  a  meridian  of  longitude  A.  is  for 
positive  longitude  (east), 

6  =  \  (23) 

for  negative  longitude  (west),  0  =  A  -f-  27r  (in  radians),  and  (?  =  A  -f  360  (in  degrees). 

The  division  of  the  northern  hemisphere  into  equal  cells,  and  the  numbering  of  the  cells, 
proceeds  as  in  the  previous  section. 

Suppose  the  northern  hemisphere  is  to  be  divided  by  parallels  of  latitude  to  make  five 
rings  M  =  5,  each  ring  with  a  width  approximately  18°.  Then  the  total  number  of  cells  I 
will  be  121  [Eq  (3)].  The  incremental  length  8r  on  the  corresponding  equal-area  map  [Eq 
(20)]  is; 


8r  =  819.08km  =  441.97nm.  (24) 

The  latitudinal  circles  that  divide  the  hemisphere  into  the  five  rings  will  correspond  [Eq 
(20)1  to  the  map  radii  r: 

r  =  m{8r) 

for  m  =  1,  3,  5,  7,  9,  11,  for  which  Eq.  (22)  gives  the  parallels: 

=  82.63°,  67.76°,  52.50°,  36.62°,  19.30°,  0(  equator) 
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For  the  location  of  Boston  <f>  =  42.5°,  A  =  —71.0°,  Eqs  (22)  and  (23)  give  r  =  2769nm, 
0  =  289®.  With  the  given  polar  map  coordinates  (r,0),  Eqs  (5)  and  (6)  will  provide  the  ring 
number  (n  +  1  =  3)  and  Eqs  (7),  (8),  (9),  and  (10)  will  provide  the  cell  number  (I  =  45). 

Problem  4:  Climatology  using  equal-area  cells. 

Suppose  it  is  appropriate  to  study  the  climate  in  the  northern  hemisphere  in  very  small 
pixels,  such  as  (5  km)^.  The  approximate  number  of  pixels  I'  will  be  given  by  the  total  area 
of  the  northern  hemisphere  divided  by  the  elenjental  area  (25  km*): 

/'=  10,187,524. 

The  number  m  of  rings  is  rounded  to  1595  [Eq.  (3)].  This,  in  turn,  modifies  the  total  number 
of  cells  to 

7  =  10,182,481 

so  that  the  area  [Eq.  (2)]  of  each  cell  6A  is  recalculated  to  be  25.01km*.  The  incremental 
length  6r  is  given  by  Eq.  (2): 

Sr  =  2.822km 

For  a  point  given  by  latitude  and  longitude,  its  cell  number  can  be  determined  using  Eqs 
(22),  (23),  (5),  (6),  (7),  (8),  and  (10).  Information  can  be  stored  in  cells  and  referenced  by 
the  cell  numbers. 

3.  CONCLUSION 

The  complexity  involved  in  dividing  a  circle  into  equal-area  pixels  is  warranted  by  the  ap¬ 
plications  of  such  a  technique.  It  can  be  applied  to  the  earth’s  spherical  surface,  and  also  to 
a  sky  dome.  The  mapping  technique  provides  a  method  for  the  equal  division  of  the  earth’s 
surface  in  preparation  for  areal  analysis  around  a  central  point,  usually  the  north  pole. 


« 
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